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Using a simulated system of temperature and concentration advected by a turbulent 
velocity, we investigated the statistical differences between active and passive scalars in 
isotropic compressible turbulence, as well as the effect of compressibility on scalar trans¬ 
ports. It shows that in the inertial range, the kinetic energy spectrum and scalar spectra 
follow the power law, and the values of the Kolmogorov and Obukhov-Corrsin 

constants are 2.06 and 0.87, respectively. The appearance of plateau in the transfer flux 
spectrum confirms the balances between dissipations and driven forcings for both ve¬ 
locity and passive scalar. The local scaling exponent computed from the second-order 
structure function exists flat regions for velocity and active scalar, while that for passive 
scalar takes first a minimum of 0.61 then a maximum of 0.73. It then shows that the 
mixed third-order structure function of velocity and passive scalar satisfies the Yaglom’s 
4/3-law. For the scaling exponent of structure function, the one of velocity and passive 
scalar mixing is between those of velocity and passive scalar, while the one of velocity and 
active scalar mixing is below those of velocity and active scalar. At large amplitudes, the 
probability distribution function (p.d.f.) of active scalar fluctuations is super-Gaussian, 
whereas that of passive scalar fluctuations is sub-Gaussian. Moreover, the p.d.f.s of the 
two scalar increments are concave and convex shapes, respectively, which exhibit strong 
intermittency at small scales, and approach Gaussian as scale increases. For the visualiza¬ 
tion of scalar field, the active scalar has ” ramp-cliff” structures, while the passive scalar 
seems to be dominated by rarefaction and compression. By employing a ” coarse-graining” 
approach, we study scalar cascades. Unlike passive scalar, the cascade of active scalar is 
mainly determined by the viscous dissipation at small scales and the pressure-dilatation 
at large scales, where the latter is substantial in the vicinity of small-scale shocklets but 
is negligible after space averages, because of the cancelations between rarefaction and 
compression regions. Finally, in the inertial range, the p.d.f.s for the subgrid-scale fluxes 
of scalars can collapse to the same distribution, revealing the scale-invariant feature for 
the statistics of active and passive scalars. 
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1. Introduction 

Scalar mixing in turbulent flows includes generation of scalar fluctuations, distortion 
of scalar interfaces, and creation of intermittent scalar gradients at small scales. It is 
of fundamental importance in a variety of fields ranging from the scattering of inter¬ 


stellar mediums in Universe (Meyer et al. 

1998 

C 

lartledge et al. 

2006 

), the dispersion 

of air pollutants in atmosphere ( 

Lu & Turco 

1994 

1995 

), to the combustion of chemi- 

cal reactions in engine ( 

Pope 1991 Warhaft 

2000 

. In the literature of fluid dynamics. 


scalar mixing is also called as ” scalar turbulence” 
was 


which has been extensively studied. It 
Obukhov (|1949|) and Corrsin (|195I|) first applied the Kolmogorov theory to passive 


scalar advected by incompressible turbulence, and derived the power-law scaling of the 
second-order scalar increment structure function 


((d.<^)2) (X (e) 


- 1 / 3 / 


. 2/3 


( 1 . 1 ) 


where Sr4> = -f r) — (j){'x.), e and are the dissipation rates of energy and scalar, 
respectively. The sign (•) denotes ensemble average. This result indicates that the cascade 
of passive scalar is similar to the classical energy cascade for velocity (Monin [l975 ). In 
particular, the scalar fluctuations are generated at large scales and transported through 
successive breakdowns into smaller scales, the cascade proceeds until fluctuations are 
dissipated at the smallest scale. By this picture, an inertial-convective range is defined 
between the largest and smallest scales. Hereafter, we call it inertial range if there is no 
ambiguity. The scalar spectrum in the inertial range is then given by 


E^ik)=Coc{e)-^^^e^)k-^/^ 


( 1 . 2 ) 


where Cqc is the Obukhov-Corrsin (OC) constant and regarded to be universal. Sreeni- 
(19961 concluded that Cqc ~ 0.68, and the experiments of passive scalar in de- 


vasan 


caying grid turbulence by Mydlarski & Warhaft (19981, for Re\ varying from 30 to 731, 
showed that Cqc = 0.45 ~ 0.55. In numerical simulation, Wang et al. (1999) studied 


the refined similarity hypothesis for passive scalar and gave that Cqc = 0.87 ± 0.10 at 


Rex = 195, and Yeung et al. (2002 2005) found that in the scalar turbulence with mean 


gradient, Cqc was around 0.67 in the range of 90 < Re\ < 700. 

Nevertheless, the Kolmogorov-Obukhov-Corrsin theory fails in describing the well- 
known anomalous scaling problem of passive scalar except for p = 3 


{\Srct>\P) (X . 


(1.3) 


It was Kraichnan (1994) first proposed a predictable model for the scaling exponent 


of passive scalar in white-noise turbulence, using a linear anstaz. Then, [Shraiman fc 
Siggia (1994) gave a series of white-noise equations for the structure functions of passive 


scalar. Although the approach of white-noise velocity is insightful and has provoked 
many theoretical and computational work, it is clear that a more realistic velocity either 


incompressible or compressible must be considered for the problem. Cao & Chen (19971 


computed passive scalar in the Navier-Stokes (NS) turbulence and developed a bivariate 


log-Poisson model. Watanabe & Gotoh (2004) studied the statistics of passive scalar in 


homogeneous turbulence and found that in the inertial range, the passive scalar scaling 
follows the bivariate log-Poisson model. 

The experimental measurements in decaying grid turbulence at high Peclet number 


(Mydlarski & Warhaft 1998) showed that the p.d.f. of passive scalar is sub-Gaussian, 


which is identical to the simulation data from Wang et al. { 1999), and Watanabe & Gotoh 


(2004). Further, the p.d.f. of passive scalar increment deviates distinctly from Gaussian 
at small scales and approaches Gaussian when scale increases. At large amplitudes, the 








































































A numerical investigation on active and passive scalars in isotropic compressible turbulence 

p.d.f. tail of passive scalar is convex, whereas that of velocity increment is concave. The 
passive scalar field has the so-called ”ramp-clifF’ structures (Chen & Kraichnan 1998 


Watanabe & Gotoh 20041, which are produced by stretching and shearing of vortices. 


Specially, turbulent mixing effectively throws out high gradient events, forming the large- 
scale ramp structures of low gradient which are reconciled with the existing small-scale 


cliff structures of high gradient (Shraiman & Siggia 20001. Thus, scalar dissipation is 


small in ramps but is large in cliffs. 

When the Schmidt number is Sc ^ 1, the sca lar fluctuations are dissipated at the 
so-called Batchelor scale rjB = Sc~^^^ri. Batchelor ( BatchelorJl959 1 predicted that in the 
viscous-convective range (i.e. ijb ^ r rj), the scalar spectrum obeys a. k~^ power law. 


which has been confirmed by many simulations such as Yeung et al. (2002 20041, Antonia 


& Orlandi (20031, and Donzis et al. (2010). On the contrary, when S'c ^ 1, the scalar 
fluctuations decay strongly at the OC scale poc = Sc~^^'^r]. Standard arguments by 
Batchelor, Howells and Townsend showed that when the Reynolds number is sufficiently 
high, the scalar spectrum follows a power law in the inertial-diffusive range (i.e. 


rj r rjoc)- Recently, this conclusion has been proved by Yeung & Sreenivasan (2012 


2014). 


The current understanding of scalar transport in compressible turbulence lags far be¬ 
hind the knowledge accumulated on the incompressible one. Previous studies involved the 
cascade of passive scalar in compressible turbulent flows were focused on the dependence 
of passive scalar trajectory on the degree of compressibility ( Chertkov et aL||19^ 1998). 
It was analyzed by Gawedzki & Vergassola (2000) that for weak compressibility, the di¬ 


rect cascade of passive scalar associating with the explosion of initially close trajectories 
takes place, which is similar to that occurs in incompressible turbulence. By contrast, 
when compressibility is strong enough, the passive scalar trajectories collapse, and there 
arise the excitation of nonintermittent inverse cascade at large scales and the suppression 


of dissipation at small scales. Celani et al. (1999) studied the dynamical role of compress¬ 


ibility in the intermittency of passive scalar for the direct cascade regime, and pointed 
out that due to the slowing down of Lagrangian trajectory separations, the reinforced 


intermittency would appear for increasing compressibility. Recently, Ni & Chen (2012) 


performed simulations of passive scalar in one-dimensional (ID) compressible NS turbu¬ 
lence, and observed that starting from certain medium scales, the transfer flux of passive 
scalar transports upscale. Nevertheless, so far there is no observation for the inverse cas¬ 
cade of passive scalar in three-dimensional (3D) compressible turbulent flows, even the 
value of turbulent Mach number reaches at 6.1 ( ]Pan fc Scannapieco|2011 ). Furthermore, 
it was found that the passive scalar scaling in supersonic turbulence computed from the 


piecewise-parabolic method (PPM) simulation accords well with the SL94 model (She 


& Leveque 1994). At low Mt, the field structure of passive scalar is similar to that in 
incompressible turbulence, while at high Mf, it is primarily dominated by large-scale 
rarefaction and compression ( Pan &: Scannapieco||20ld ). 

Gontrasting to passive scalar, an active scalar always provides an impact on velocity 
by way of an extra term, i.e. buoyancy, added to the velocity equation, resulting in a 
two-way coupling between the velocity and scalar fields. Thus, it belongs to nonlinear 
problems. Celani et al. ( 2002a|b 2004) studied the university and scaling of active and 
passive scalars by simulating four different numerical systems, and obtained the results 
that when the statistical correlation between active scalar forcing and advecting velocity 
is sufficiently weak, the active and passive scalars share the same statistics; however, 
when the statistical correlation is strong, the two scalars behave quite different. In the 


context of shell model, Ching et al. (2003) argued that the different scaling behaviors 


are attributed to the fact of the active scalar equation possesses additional conservation 
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laws, while in the passive scalar equation the zero mode problem is a secondary factor. 


In the thermal convection experiment, Zhou & Xia (2002) showed that the active and 


passive scalars share some similar statistical properties, such as the saturation of scaling 
exponent and the log-normal distribution of front width. 

In this paper, we regard the temperature advected by a compressible turbulent velocity 
as an active scalar. Without an additional term to the velocity equation, the tempera¬ 
ture brings substantially nonlinear feedback to the velocity through the way of density 
fluctuations, which fits the general spirit of active scalar. We carry out a numerical in¬ 
vestigation on the turbulent temperature and concentration transports. Our attention is 
focused on the statistical differences between active and passive scalars, as well as the 
effect of compressibility on scalar fields. The system is simulated on a 512^ grid using 
a novel hybrid method. This method utilizes a seventh-order weighted essentially non¬ 


oscillatory (WENO) scheme ( 

Balsara & Shu 

to 

0 

0 

0 

in shock region and an eighth-order 

compact central finite difference (CCFD) scheme 

Tele 

19921 in smooth region outside 


shock. In addition, a new numerical hyperviscosity formulation is proposed to improve the 
numerical instability in simulation without compromising the accuracy of the method. As 
a result, the hybrid method has spectral-like spatial resolution in computing highly com¬ 
pressible turbulence. In comparison, though the PPM method can simulate compressible 
flows at higher turbulent Mach number, the absence of viscous dissipation brings a series 
of unphysical problems. The details of the hybrid method has been described in |Wang| 
jet al. (2010). We hope that our comprehensive study will advance one’s understanding 


of the statistics of active and passive scalars in the small-scale compressible turbulence. 

The remainder of this paper is organized as follows. The governing equations and 
basic parameters, along with the computational method, are described in Section 2. The 
fundamental statistics of the simulated system is provided in Section 3. The probability 
distribution function and structure function scaling are analyzed in Section 4. In Section 
5, we describe the statistical behaviors of field structures. The discussions on the cascades 
of active and passive scalars are given in Section 6 . Finally, we give the summary and 
conclusions in Section 7. 


2. Governing equations, simulation parameters and methods 

We consider a dynamic system of temperature and concentration transports by a tur¬ 
bulent velocity field under the idea gas state. The system is driven and maintained 
in a stationary state by the large-scale velocity and passive scalar forcings. Similar to 


Samtaney et al. (2001), we shall introduce a set of reference scales to normalize the 


hydrodynamic and thermodynamic variables. Since these variables together contain five 
elemental dimensions, we first introduce five basic reference scales as follows: L for length, 
Po for density, U for velocity, Tcq for temperature, and Coq for concentration. Then we 
obtain the reference sound speed cq = y/jRTeo, kinetic energy per unit volume poU^/2, 
pressure poRTeo, and Mach number M = U/co, where 7 = Cp/Cy is the specific heat 
ratio, R = Cp — Cy is the specific gas constant, and Cp and Cy are the specific heats at 
constant pressure and volume, respectively. By adding a reference dynamical viscosity po, 
thermal conductivity Kq, and molecular diffusivity xo^ derive three additional basic 
parameters including the reference Reynolds number Re = poUL/p^, Prandtl number 
Pr = poCp/kq, and Schmidt number Sc = I'ojxo. In our simulation, the values of 7 , 
Pr and Sc are set as 1.4, 0.7 and 1.0, respectively. Then the flow system is completely 
determined by the two parameters of M and Re. 

Based on the above reference scales, the dimensionless form of the governing equations 
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for the simulated system are written as 


d d \ d 

717(^“0 +p5,,hM^] = 


dt 


Re dxj 


( 2 . 1 ) 

( 2 . 2 ) 


f) f) 1/9 f)'T' p 1/9 

[(£ +p/7^")m,] = - a + pJ^jUj, ( 2 . 3 ) 


a dxj ^ dxj Re dxj 
d I ^ s d r, „ -1 1^ / dCo. 

Si + 37 = ;s57<™ &7> + 

p = pTe, 


(2.4) 

(2.5) 


where a = PrRe{j— 1)M^ and f3 = ScReipf— 1 ) 7 . The primary variables are the density 
p, velocity Ui, pressure p, temperature Te, and concentration Co. Rj = J2m=i ^i^m) x 
exp(ikmx) +C.C. is the large-scale velocity forcing, where is the Fourier amplitude and 
is perpendicular to the wavenumber vector ( Ni &: Chen|2012 Ni et al. 20131. Similarly, the 
large-scale passive scalar forcing is 5 = ^^(kni) exp(ikmx ) -|- c.c., and the F ourier 

amplitude S is also perpendicular to the wavenumber vector (Ni & Chen 2012 1 . The 


inclusion of the large-scale cooling function A in the energy equation is to remove the 
accumulated internal energy at small scales, and its detail can be found in |Wang et al. 
(2010). The viscous stress aij and total energy per unit volume £ are defined by 


, dui dui 




£ = 




( 2 . 6 ) 


(2.7) 


(7 — 1 ) 7 M^ 

where 6 = duk/dxk is the dilatation, a flow variable that measures the local rate of 
expansion or compression. Here the dimensionless dynamical viscosity /r, thermal con¬ 
ductivity K and molecular diffusivity x are considered to vary with temperature through 
the Sutherland’s law ( Sutherland||1992 ) as follows 

1.4042Tei'5 


P,K,X = 


Te + 0.4042 ’ 


( 2 . 8 ) 


The system is solved numerically in a cubic box with periodic boundary conditions. 


by adopting a new computational approach (Wang et al. 2010). To obtain statistical 
averages of interested variables, a total of eighteen flow fields spanning the time period 
of 6.15 < tjx ^ 10.51 were used, where r = 1.17 is the large-eddy turnover time, and the 
system achieves the stationary state at the time of t > At. Instead of the reference Mach 
number M and Reynolds number Re, in simulation the flow is directly governed by the 


turbulent Mach number Mt and Taylor microscale Reynolds number Rex (Samtaney et 

^[200T1 ) 


Mt = M— 

(Vre)’ 


Re\ = Re 


u'Hp) 

V3(m)’ 


(2.9) 


( 2 . 10 ) 


where u' = •\/(uf~+u|~+u|) is the r.m.s. magnitude of velocity, and A is the Taylor 
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Table 1. Simulation parameters and resulting flow statistics. 


Grid 

Mt 

Rex 

M' 

kmaxf] 

Lf 

Lf/V 

A 

T S 3 

e’ 

512® 

1.02 

178 

1.04 

2.79 

1.50 

135 

0.27 

1.17 

-1.7 

5.8 

u' 

e'/x' 

u' 


Ucju's 

P' 

(^) 


{ec)/{e} 

Ek 

El 

15.6 

0.37 

2.21 

2.16 0.50 

0.23 

0.29 

0.56 

86.8% 

16.0% 

2.31 

8.82 


microscale, which is defined as 


,/2 


A = 


{{dui/dxi)'^ + {du2/dx2Y + [du^/dx^Y)' 


( 2 . 11 ) 


In our simulation, by setting M = 0.45 and Re = 500, we obtained the time-average 
values of Mt = 1.02 and Re\ = 178. Although there is hardly a real instance with 
such values of Mt and Re\, the investigation reveals certain universal features of inertial 
ranges for compressible flows at different Reynolds numbers. 


3. Fundamental statistics of simulated system 

Table lists some overall statistics of the velocity field. The resolution parameter is 
kmaxV = 2.79, where kmax and i] are the largest solved wavenumber and the Kolmogorov 
scale T] = [{fj,/{Rep)Y/ < ejp and their values are 256 and 0.011, respectively. The 

integral scale T/ is computed by 


= f Y^dk. 


2u'2 


k 


(3.1) 


where E{k) is the kinetic energy spectrum per unit mass, namely. 


OO 

J E{k)dk = u'^/2. (3.2) 

0 


The ratio Lf/rj is 135, representing the range of scales in the velocity field. The r.m.s. 
magnitudes of dilatation and vorticity are computed by 9' = dj ( 0 ^), ix' = dj+x\-\- w|), 
with respective values of 5.8 and 15.6. Then the ratio 9'/x' is 0.37, indicating that the 
compressibile effect makes a significant contribution to the statistics of velocity gradient 
tensor. This character is further demonstrated by the ensemble average of the velocity 
derivative skewness, which is defined by 


^ + + (i5)T/3 

[«fe)" + (5g)“ + (iS)T3l“'" 


The value of S '3 is equal to —1.7, much greater than typical values of —0.4 ~ —0.6 
observed in incompressible turbulence (Ishihara et al. 2007). This deviation is mainly 
caused by the presence of the small-scale shocklets in compressible turbulence. 

As for the scalar fields, their overall statistics are summarized in TableThe effective 
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Table 2. Simulation parameters and statistics of scalars. 


Scalar 

kmaxlJOC 


Lj^/goc 


Ssfj) 


(«</>) 


Te 

3.64 

0.79 

54 

0.16 

0.13 

0.12 

0.06 

0.51 

Co 

2.79 

0.96 

86 

0.23 

-0.11 

2.25 

1.45 

2.48 


resolution parameters are kmaxVOC = 3.64 for temperature and 2.79 for concentration. 
The integral scale for scalar is computed by 

OO 

= <«) 

0 


Hereafter (p is denoted as Te or Co whenever there is no special illustration. E^{k) is the 
scalar spectrum per unit mass, and its definition is similar to Equation (3.2). is the 
r.m.s. magnitude of scalar. The ratio Lf^/rjoc for temperature and concentration are 
54 and 86, respectively. It shows that the ranges of scales for the two scalars are smaller 
than that for the velocity. Furthermore, the ensemble average of the skewness of scalar 
derivative is defined by 


^30 


+ & + (&W3 
[«&)" +(ft)"+ (g)")/3]’'" 


(3.5) 


The values of Ss^i) are 0.13 for temperature and —0.1 for concentration, significantly 
smaller than S 3 in magnitude. Besides, the Taylor microscale X^, defined similarly to A, 
are 0.16 for temperature and 0.23 for concentration. 


The application of Helmhotlz decomposition (Erlebacher & Sarkar 1993) to velocity 
field gives that 

Ui = + Ui,c, (3.6) 

where Ui^s is the solenoidal component and satisfies dui^s/dxi = 0, and Ui^c is the com¬ 
pressive component and satisfies £ijkduk,c/dxj = 0, and Sijk is the Levi-Civita symbol. In 
Table we compile the relevant statistics associated with the solenoidal and compressive 
processes. The r.m.s. magnitudes of solenoidal and compressive velocity are = 2.16 
and = 0.50. This leads the ratio of to be 0.23. We also note that 9'/uj' is larger 
than u'^ju'g, indicating that the compressible effect plays a more significant role at the 
small scales. The r.m.s. magnitudes of scalars presented in Table are Te' = 0.12 and 
Co' = 2.25. 

Let us now focus the governing equation of the kinetic energy pUiUi/2. According to 
Andreopoulos et al. (20001, it can be written as follows 


D 


/I \ 

d / 

1 \ 


dxj V 



1 A* 


4 M .2 c, duj 2\ , X- /Q '7\ 


where D/Dt = d / dt+Ujd / dxj is the material derivative. On the right-hand side (R.H.S.) 
of Equation (3.7), the viscous dissipation rate e = aijSij/Re is divided into three parts: 
the solenoidal dissipation rate Cg = {p,/Re)u}iUJi, the dilatation dissipation rate Cc = 
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Figure 1. Compensated spectra of kinetic energy and its two components, and compensated 

spectra of scalars. 


{4:/3){fj,/Re)9^, and the mixed dissipation rate = {2^/Re)[{dui/dxj){duj/dxi) — 0^], 
which represents the contribution to dissipation rate from the non-homogeneous compo¬ 
nent of compressible flow, and is zero in the incompressible limit. We find that (es)/(e) = 
86 . 8 % and (ec)/(e) = 16.0%, meaning that the mixed part gives a negative contribution 
to the viscous dissipation rate. The ratio Ek/Ej is 0.26, close to the estimated value of 
Ek/Ej « (7 — l)^Mf /2 = 0.29, where Ek = (pUiUi) j2 and Ej = (p )/[(7 — 1 ) 7 M^]. In 
addition, the ensemble-average value of the injected energy is (pEiUi) = 0.55, implying 
that there is a main balance between the velocity forcing and viscous dissipation. As 
long as the velocity and temperature helds are statistically stationary, one must have 
{pFrUi) = (A). 

By Equations (2.2)-(2.3), we derive the governing equation of temperature as follows 


Wt 


(pTe) 


^}ipTe)uj] = {-f-l)p9 + 


7 d . dTe 
Pe dxj ^ dxj 


a"f 

PeRe^ 


duj 

dxj 


a7 


A, (3.8) 


where Pe = PrRe is the Peclet number. The dissipation rates of temperature and con¬ 
centration are defined by 

cte = KidTeldx^)^, (3.9) 

e™ = x{dColdxjf. (3.10) 

It shows in Table that (cte) and (cco) are 0.06 and 1.45, respectively, and {CoS) is 
1.44. This confirms the approximative balance between the passive scalar forcing and 
molecular dissipation. The values of the temperature and concentration variances per 
unit volume are Et^ = 0.51 and E^o = 2.48, where Ef^ = (pTe^)/2 and E^o = (pC'o^)/2. 

Figure]^ shows the compensated spectra of the kinetic energy E{k) and its two com¬ 
ponents Es{k) and Ec{k), as well as those of the scalars Ete{k) and Eco{k). The inertial 
ranges of these spectra are clearly identified in the range of 0.05 < fciy < 0.18. We observe 
that the Kolmogorov constant of E{k) and OC constants of Ete{k) and Eco{k) read from 
the spectra are 


Cic = 2.06 ± 0.03, Coc = 0.08 ±0.01, Coe = 0.87 ± 0.02. (3-11) 


Here Ck is slightly higher than the typical values observed in incompressible turbulent 
flows (Wang et a7||1996), and Cgo is in well agreement with that from passive scalar in 
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Figure 2. Normalized flux spectra of kinetic energy and scalars. 


incompressible turbulence provided by Wang et al. (1999). However, Cq(j is roughly one 
order of magnitude smaller than Cqq. The reason is due to the fact that the transfer flux 
of active scalar is more tightly associated with the compressive component of velocity, 
which contributes a small fraction to the kinetic energy. The spectra of kinetic energy 
and its solenoidal component almost overlap, except at high wavenumbers where the 
compressive component dominates the energy content. Furthermore, for the compressive 
component, the compensated energy level in the inertial range is roughly one order of 
magnitude smaller, but decays more slowly in the dissipative range. For scalars, the con¬ 
centration spectrum is close to the kinetic energy spectrum and its solenoidal component, 
while the temperature spectrum resembles the compressive component. 

Similar to Watanabe & Gotoh (2004 2007), we define the transfer functions of the 


kinetic energy and scalars in wavenumber space as follows 

{^^+2^lk^)E{k) = T{k)+S{k), 

E(j,[k) = T^(k) -\- ^,f,{k), 


(3.12) 

(3.13) 


where S{k) and ^^{k) are the terms of the velocity and passive scalar forcings in wavenum¬ 
ber space. Then the transfer fluxes of the kinetic energy and scalars across the wavenum¬ 
ber /c are 


OO OO 

n(A:) = J T{k)dk, H^)/;) = J T^{k)dk. 


(3.14) 


k k 

In Figure we plot the spectra of transfer fluxes for the kinetic energy and scalars 
against the dimensionless wavenumber kp, normalizing by their ensemble averages of 
dissipation rates. It is found that n(fc) and nco(fc) are close to unity in the ranges of 
0.02 ^ kr] ^ 0.04 and 0.02 ^ kp ^ 0.05, respectively, over which E{k) and E^oik) are 
in the inertial range and thus have the k~^^^ power law. By contrast, in the range of 
0.05 ^ krj ^ 0.09, Ilteik) is only 24% of {cte)- This indicates that the dominance of the 
rarefaction and compression motions leads the flux transfer of temperature to being not 
follow the Kolmogorov picture. The shift towards higher wavenumbers is due to that the 
OC scale of temperature is larger than the Kolmogorov scale. Furthermore, the constancy 















10 


Qionglin Ni 



Figure 3. The p.d.f. of the ratio of the normalized dilatation to vorticity magnitude. Inset: 
the p.d.f. of the normalized dilatation. 



of the transfer-flux spectra in certain wavenumber ranges indicate that the velocity and 
scalar fields are really in the equilibrium states. 

The p.d.f. of the ratio of the normalized dilatation to vorticity magnitude is shown 
in Figure]^ The p.d.f. is negatively skewed, exhibiting that the events of large negative 
ratio is substantial, even though 9'/to' has a small value of 0.37. Similar to those of 
dilatation, the p.d.f. tails of the ratio are concave. In the inset we plot the p.d.f. of the 
normalized dilatation with strongly negative skewness. This feature has already been 


observed in weakly and moderately compressible turbulent flows (Porter et al. 2002 


Pirozzoli & Grasso 20041. It was claimed by Pirozzoli fc Gras's^ (2004) that the p.d.f. 


becomes more skewed towards the negative side as the turbulent Mach number increases. 
This point has been confirmed in our simulation. In Figure we plot the p.d.f. of the 
local Mach number Mioc- Glearly, a substantial portion of the velocity field is supersonic, 
where the maximum of Mioc is 3.32. The p.d.f. peaks at Mioc = 0.86. It reveals that 
in compressible turbulence, the small-scale shocklets generated by velocity fluctuations 
distribute randomly at different scales. 

The structure functions of the longitudinal velocity increment 5rU = u(x -I- r) — u(x), 
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Figure 5. Second-order structure functions of velocity and scalars, and mixed second-order 
structure functions of velocity-scalars, as functions of r/g. Inset: local slopes as functions of 

r/g. 


and the scalar increment Srfj) = f) ~ defined by 


SuA'^) = (|^rM|P), S'0,p(r) = 


(3.15) 


where r is the separation distance between two points and p is the order number. The 
Kolmogorov theory predicts that the structure functions scale as In Figure we 
plot the second-order structure functions of 5rU, SrTe and SrCo against the normalized 
separation distance r/g. It shows that for SuA’’’) 5'te,2(u), the regions corresponding 

to the flat local scaling exponents are observed in the range of 50 ^ r/g ^ 110. By 
contrast, for ^co, 2(^)1 in the same range there are two different regions, in which the local 
scaling exponent takes first a minimum and then a maximum, and crossover occurs in 
the range of 65 ^ r/r; ^ 95. Further, we plot the mixed second-order structure functions 
of velocity and scalars, with the relevant definition as follows 

(3.16) 


It shows that the behaviors of SmteA'^) ^nd SmcoA"^) similar to those of SuA'^) 

< 5 ( 6,2 (u), namely, for the local scaling exponents, there are neither local minima nor local 
maxima. Thus, we speculate that throughout the range of 50 < r/g < 110, the transfer 
functions of both temperature and concentration at a given order number have single 
scaling exponents. 

In the inset we plot the local scaling exponents as functions of r/g, defining as 


ZuAf) 


d\og[Su,2{r)] 
d\og(r/g) ’ 


^ 0 , 2(0 


dlog[50,2(r)] 

dlog(r/? 7 ) 


(3.17) 


. = dlog[5„^0.2(r)] 


(3.18) 


dlog{r/g) 

In the range of 60 ^ r /77 ^ IIO, Zu, 2 ii") and z^teA"^) almost constant and have the 
around values of 0.67 and 0.69, respectively. For Zte, 2 (c) and Zmco, 2 ir), their constant 
values of 0.77 and 0.72 are found in the ranges of 90 ^ r/g ^ 110 and 60 ^ r/g ^ 90, 
respectively. In terms of Zco, 2 {i’), it first decreases from small scales, and achieves a local 
minimum of z/!//^ — 6-61 at r/g = 55, then quickly increases and reaches a local maximum 


of , 


= 0.73 at r/g = 100, after that, decreases again. To some extent, the behaviors 
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of the second-order structure functions and their local scaling exponents in this study are 
similar to those observed in the ID compressible turbulence ( Ni &: Chen|2012 ). However, 
so far there are still unclear in two aspects: ( 1 ) whether a flat region in Zco, 2 {j’) appears 
at the level of the local minimum or maximum as the Reynolds number increases; and 
( 2 ) whether the local minimum and maximum features of 200,2 (?') disappears when the 
Schmidt number 5'c 1. 

In the final of this section, we focus our attention on the applicability of the Kol¬ 
mogorov’s 4/5 and Yaglom’s 4/3-laws to compressible turbulence. The Kolmogorov’s 
4/5-law is an asymptotically exact result for incompressible turbulent flows, and can be 
obtained from the Karman-Howarth equation ( Pope|[2000 ). When the Reynolds number 
is sufficiently high, the velocity field of a stationary incompressible turbulence obeys the 
following equation 


{{Sruf) = + Gv-^iiSruf) + f. (3.19) 

In the range oir] r Lf, the second and third terms on the R.H.S. of Equation (3.19) 
vanish, and then it yields to the 4/5-law 


((^ru)^) = -^(e)r. 


(3.20) 


Similarly, the Yaglom’s 4/3-law is the asymptotically exact result for incompressible 
scalar turbulence, and was derived by Yaglom in 1949 ( Yaglom||l949| . It reads as 

d 


{5rU{6r4>) ) = + 2v^ — {{5r4>) ) +/ 0 , 


(3.21) 


where denotes the molecular diffusivity. When rjB ^ r Equation (3.21) is 

reduced to the 4/3-law 

{5ru{5r(l)f) = - (3.22) 

It should be pointed out that the smallest scale for the 4/3-law to hold is 775 rather 
than 77 . In this study the Batchelor scales for temperature and concentration are rjB = 
Sc~^l'^rj = 77 and r]B = Pr~^/'^ri « I. 277 , respectively. 

In Figurej^we plot the third-order velocity structure function —Su^sir) and the mixed 
third-order velocity-scalar structure functions —Smte,?,{'>') and —Smco, 3 ir), against the 
normalized separation distance r/rj. For —Su, 3 {r) and —Smco, 3 {‘r), there are plateaus in 
the ranges of 18 ^ r /77 ^ 35 and 60 ^ r /77 < 105, respectively. Contrarily, there is 
no clear plateau observed in —Smte, 3 ir). Instead, it peaks at the around value of 0.18, 
sharing a similar explanation with the spectrum of temperature transfer flux. The inset 
shows the details of flat regions for 10 ^ r /77 ^ 140. The level values computed from 
the plateaus are 0.78 and 1.39, which are very close to 4/5 and 4/3, respectively. This 
concludes that the velocity and passive scalar in our simulation satisfy the exact relations 
of structure functions derived from incompressible flows. 


4. Probability distribution function and high-order statistics 

As we known, the intermittency in turbulence can be characterized by the behaviors 
of the p.d.f. tails of velocity and scalar increments. Figure shows the time-average 
one-point p.d.f.s of the normalized fluctuation components of velocity, temperature and 
concentration, where u/ Te' and Co' have been referred in Section 3. The p.d.f. with 
Gaussian distribution is plotted for comparison. We find that at small amplitudes, the 
p.d.f.s of velocity and temperature are very close to Gaussian, while that of concentration 
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Figure 6 . Compensated third-order structure function and mixed third-order structure func¬ 
tions of velocity-scalars, as functions of r/g, respectively. —Su,3- circles; —Smte,3'- squares; 
— Smco,3'- diamonds. Inset: the same plot of —Su,3 (solid line) and —Smco,3 (dotted line) in 
the range of r/g = 10 ~ 140. 



Figure 7. The p.d.f.s of the normalized fluctuations of velocity and scalars. The dash-dotted 

line represents Gaussian. 


occurs small oscillations. The reason is that in turbulence the fluctuations of concentra¬ 
tion are usually more intense than those of velocity and temperature. At large amplitudes, 
the p.d.f. of temperature decays smoothly and more slowly than Gaussian, and thus is 
called as super-Gaussian. This feature is also found from the scalar in 3D Kraichnan- 
model flows (Balkovsky fc Lebedev|19^ Mydlarski fc Warhaft|1998 Shraiman & Siggia 


2000 Gelani et a/.||2002a|b |2004[ ). By contrast, the p.d f.s of velocity and concentration 

are called as sub-Gaussian, because they decay smoothly and more quickly than Gaus¬ 


sian, similar to the passive scalar in incompressible turbulence (|Mydlarski & Warhaft 
1998 Watanabe fc Gotoh||2004 1. 


The time-average one-point p.d.f.s of the normalized longitudinal gradients of velocity, 
temperature and concentration are shown in Figure]^ where ^ = dujdx^ ( = dTe/dx 
and = dCo/dx. It shows that the three p.d.f.s all deviate drastically from Gaussian 
(dash-dotted line), exhibiting strong intermittency. The p.d.f. of ^ is negatively skewed. 
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Figure 8 . The p.d.f.s of the normalized longitudinal velocity and scalar gradients. Inset: log-log 
plot of the p.d.f.s of the negative-value gradients, where the numerals represent the values of 
the line slopes. 


while those of C and are approximately symmetric. In the negative gradient regions, 
the p.d.f. tails of ( and c are longest and shortest, implying strongest and weakest in 
intermittency, respectively. However, in the positive gradient regions, the relative inten¬ 
sity of intermittency between ^ and c reverses. Furthermore, in the vicinity of varnishing 
gradient, the p.d.f.s of C and c peak higher than that of <^. In the inset we present the 
log-log plot of the large negative gradients versus their p.d.f.s. It shows that for C 
and <j, the values of the power-law exponents are —3.8, —3.3 and —3.5, respectively. The 
major contributions of these p.d.f. tails are caused by small-scale shocklets rather than 
large-scale shock waves. Here we notice that the exponent value of —3.8 for the p.d.f. of ^ 
is a bit larger than that for the p.d.f. of velocity gradient, i.e. —3.5, in Burgers turbulence 
( E fc EijndenlpoM 2000 Bee et ar]|2000 ), but is substantially larger than that for the 
p.d.f. of dilatation, i.e. —2.5, in compressible turbulence ( Wang et a/.||2012a ). 

In Figurej^we plot the p.d.f.s of the normalized longitudinal velocity increment Srujasu 
and its components SrUs/asu, and SrUcjcrsu^, at the normalized separation distances of 
rfrj = 1, 4, 16, 64 and 128, where asu, and asu^ are the standard deviations of SrU, 
SrUs and SrUc, respectively. We find that in Figure a) at small scales (i.e. r/rj = 1), 
the left p.d.f. tail of 6rU is much longer than the right one, indicating strongly negative 
skewness. This feature is similar to the p.d.f. of longitudinal velocity gradient shown in 
Figure]^ The p.d.f. of SrU gets more narrow and symmetric as r/ry increases, and it 
recovers to Gaussian without intermittency at r /77 = 128. At large amplitudes, the p.d.f. 
tails are concave, which is in agreement with those observed in the experiments and 
simulations of incompressible turbulence. As for the solenoidal component 6rUs shown in 
Figurej^b), at small scales, on each side, the p.d.f. tails are slightly negative asymmetry 
and longer than Gaussian. When r/rj increases, they quickly become Gaussian. We notice 
that at r /77 = 1, the p.d.f. of SrUg is close to that of SrU provided by Ishihara et al. (2009) 
(circles), implying that the dynamics of solenoidal-component velocity is similar to that of 
velocity in incompressible turbulent flows. On the other side, the p.d.f. of the compressive 
component 5rUc shown in Figure |^c) displays strong intermittency at small scales, and 
thus, is like that of 5rU in Burgers turbulence (deltas). As r /77 increases, it approaches 
Gaussian as well. Besides, at large amplitudes the p.d.f. tails of the velocity component 
increments are also in concave shape. 
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Figure 9. The p.d.f.s of the normalized longitudinal velocity increment and its two components 
at various r/g values. The black solid line represents Gaussian. The symbols of circles and deltas 
in (b) and (c) are for the p.d.f.s of the normalized longitudinal velocity increment at r/g = 1 in 
incompressible turbulence and Burgers turbulence, respectively. 



Figure 10. The p.d.f.s of the normalized scalar increments at various r/g values. The black 

solid line represents Gaussian. 


The p.d.f.s of the normalized scalar increments of SrTejasTe and 5rCo/asco are shown 
in Figure where cr^Te and asco are the standard deviations of 6rTe and 6rCo, re¬ 
spectively. Obviously, these p.d.f.s are almost symmetric, even at small scales such as 
r/g = 1. When r/g increases, they get narrower and approach Gaussian. It shows that 
on both sides the p.d.f. tails of S^Te are significantly longer than those of SrCo, indi- 
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Figure 11. Skewness of velocity and scalar increments as functions of r/rj. 


eating that in intermittency SrTe is much more intense than SrCo . Furthermore, at 
large amplitudes, the p.d.f. tails of temperature increment are concave, similar to those 
of longitudinal velocity increment, while those of concentration increment are convex. 
This feature has been found in the previous simulations of passive scalar advected by 
incompressible turbulence ( Watanabe &: Got^|2004 ). 

The change of the distributions of velocity and scalar increments can be characterized 
by their skewness and flatness as well. Those for the longitudinal velocity increment and 
scalar increment are defined as follows (Watanabe fc Gotoh||2004| 




{(Sruf) 

(((5,u)2)3/2’ 


Ku,A(r) = 


{{SrUW 


= 


K4,A^) = 


(4.1) 


(4.2) 


((< 5 ,</.) 2 ) 3 / 2 ’ (( 5 , 0 ) 2 ) 2 - 

In Figure [TT] we plot the skewness of 6rU, 6rTe and SrCo against the normalized separation 
distance r/rj. It shows that in a wide scale range, KuAr) ^'^d Kco, 3 ir} are negative and 
their local minima are around — 2.1 and — 0 . 1 , respectively, whereas KteA'^) i® positive 
and its local maximum is around 1.4. As r/rj increases, the magnitudes of skewness fall 
and approach zero. We find that the Gaussianities of 5rU, 6rTe and SrCo begin to appear 
even the scale is not very large. Throughout scale ranges, the magnitude of Kco, 3 {r) is 
always smaller than that of Ku^sir), implying a weaker degree of Gaussian departure. 

The flatness of S^-u, SrTe and SrCo against the normalized separation distance r /77 is 
plotted in Figure 12 It shows that when r/rj decreases, the three flatness all increase 
from the Gaussian value of 3.0, meaning that the intermittency of the flow increments 
emerge and enhance. In a wide scale range, KuA’’^) i® smallest, and thus, is weakest in 
intermittency. Compared to KcoAAj T^teAA i® Is^rger at small and moderate scales. 
However, the relation reverses in the range of 20 ^ r/rj ^ 90, indicating that the in¬ 
termittency of SrCo surpasses that of SrTe. This behavior has not been announced in 
the text of probability distribution function. The detailed variations of flatness in the 
range of 100 ^ r/rj ^ 400 are shown in the inset. At sufficiently large scales, KuA^) 
and KcoA’’^) below 3.0 and thus are sub-Gaussian, while KteA"^) i® above 3.0 and 
thus are super-Gaussian. This feature is consistent with the behaviors of the p.d.f.s of 
fluctuations at large amplitudes shown in Figure 
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Figure 12. Flatness of velocity and scalar increments as functions of r/g. Inset: the same plot 

in the range of r/g — 100 ~ 400. 



Figure 13. Scaling exponents as functions of the order number p. Velocity: solid circles; tem¬ 
perature: solid squares; concentration: solid diamonds; velocity-temperature: open squares; ve¬ 
locity-concentration: open diamonds. The solid deltas are for the passive scalar scaling given 
by Watanabe & Gotoh. The dashed and dash-dotted lines are for the K41 and SL94 scalings, 
respectively. 


In Figure [T^ we plot the scaling exponents of the structure functions of velocity and 
scalar increments as well as the mixed structure functions of velocity-scalar increments, 
as functions of the order number p. These scalings are computed by taking averages on 
the local scaling exponents, where the curves are in the flat regions with finite widths. It 
shows that at large order numbers, Zu,p > Zco,p > zte,p, indicating that the velocity and 
temperature fields are weakest and strongest in intermittency, respectively. As for the 
mixed scaling, we obverse that Zmco,p locates between z^^p and Zco,p, similar to that in 
incompressible turbulence. Nevertheless, it is interesting to find that Zmte,p locates below 
Zu,p and Zte,p- This phenomenon was previously appeared in 3D weakly compressible 


turbulence (Benzi et al. 20081 and ID compressible turbulence (Ni & Chen 20121. For 


comparison, we plot the velocity scaling given by the K41 ( Kolmogorov||I94T ) and SL94 


(She & Leveque 19941 models, and the passive scalar scaling provided by Watanabe 
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Figure 14. Mixed scaling exponents and their components as functions of the order number p. 
The solid symbols of circles, squares and diamonds are for velocity-temperature and its compo¬ 
nents, while the open symbols of circles, squares and diamonds are for velocity-concentration 
and its components. The symbols of deltas and gradients separately represent the scaling ex¬ 
ponents of the mixed velocit y-temperature an d velocity-concentration structure functions from 
ID compressible turbulence (Ni & Chen|2012l. 


& Gotoh (20041. Obviously, Zu,p is close to the SL94 model and Zco,p is smaller than 


Watanabe’s result. In our simulation, there is no saturation observed for the scalings of 
temperature and mixed velocity-temperature. Here we point out that the above results 
are obtained by using long-time averages (eighteen large-eddy turnover time periods). In 
fact, for each single time frame, the oscillation intensity of passive scalar scaling is much 
larger than those of passive scalar and velocity scalings. 

The Helmholtz decomposition on the mixed velocity-scalar scalings are depicted in 
Figure 14 It shows that Zmco,p and its two components Zmcos,p and Zmcoc,p are almost 
identical, and locate far away from Zmco,id- By contrast, Zmte,p is a bit smaller than its 
solenoidal component Zmtes,p, and the compressive component Zmtec,p nearly collapses 
onto Zmte,id- This indicates that the scaling for the compressive component of the mixed 
velocity-temperature is similar to the Burgers scaling, which saturates at order number 
p ^ 5. 


5. Behaviors of scalar field structures 

So far we have focused on quantifying the fundamental statistics, probability distribu¬ 
tion function and high-order statistics. In this section, we shall explore the flow structures. 
We hope that the structural analysis will complement the statistical results and provide 
a better understanding of the active and passive scalars in compressible turbulence. 

5.1. Field structure 

In Figurej^we plot the two-dimensional slices (x-y plane) through the active and passive 
scalar fields simultaneously in the position oi z = tt/2. The temperature field is depicted 
in the left panel, it shows that the small-scale ’’cliff’ structures are separated by the 
large-scale ’’ramp” structures. The temperature fluctuations are small in the divided 
structure regions, however, they quickly increase when across the interfaces between 
cliffs and ramps. To some extent, this picture is similar to passive scalar in incompressible 
turbulence producing by stretching (Shraiman & Siggia 20001, which rearranges scalar 
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Figure 15. Two-dimensional slices at ^ = 7r/2. (a) Temperature; (b) Concentration. 


gradients and brings fluid elements with very different scalar level next to each other. By 
contrast, the concentration field shown in the right panel displays to be quite different. 
The concentration fluctuations seem to be dominated by the motions of rarefaction and 
compression, causing by shock fronts. Compared with the PPM simulations provided 
by Pan & Scannapieco (2010), the concentration field in this study is similar to that 
in the = 6.1 flow, showing the effect of shock waves. Nevertheless, we will show in 
the following section that the cascade of passive scalar is indeed mainly governed by the 
stretching of turbulence. 

In Figure the 3D isosurfaces of the magnitudes of dilatation 9 = —39', vorticity 
oj = 1.8a;', temperature Te = l.lTe', concentration Co = —OACo', l.&Co' are displayed 
in a same 128^ subdomain (i.e. 1/64 of the full domain). These isosurfaces are colored 
based on the local r.m.s. magnitudes, where the yellow isosurfaces stand for the large-scale 
shock wave with a thickness typically much smaller than the inertial scales, surrounding 
by the small-scale shocklets. In Figure [T^ a), the green isosurfaces are the vortices under 
random distribution. When moving away from the shock wave, the geometric shape of 
vortices changes from sheetlike to tubelike. In Figure[l6j)b), the green sheetlike isosurfaces 
are temperature, exhibiting as wrinkled curves around shock wave. We notice that the 
surfaces of dilatation and temperature intersect with each other in certain positions. The 
isosurfaces of the positive and negative r.m.s. magnitudes of concentration are depicted 
in Figure [T6 |)c) by the colors of blue and red, respectively. It shows that the concentration 
surfaces are sheetlike as well, and the curve surfaces of oppositive sign basically appears 
on each side of the shock front. 

The p.d.f.s and conditional p.d.f.s of cosines of angles between vorticity and scalar 
gradients are shown in Figure [T7| Unlike the positive alignments between vorticity and 
vortex stretching vector in turbulent flows (Kholmyansky et al. ]MT||Wang et aZ.|20I2bD , 
the scalar gradients prefer to be perpendicular to the vorticity. In other words, the p.d.f.s 
get their maximal values in the vicinity of the tt/2 angle. For each dilatation range, the 
p.d.f. and conditional p.d.f.s of concentration are much steeper than those of temperature, 
implying that the surfaces of passive scalar and vorticity are probably tangent to each 
other. In addition, we find that the p.d.f.s of the two scalars conditioned at each diltation 
level have similar relative relations. 


5.2. Dissipation field structure 

The normalized dissipation spectra of kinetic energy and scalars are plotted in Figure [TSl 
The maximum dissipation spectrum of kinetic energy is about 2.18 and occurs at r/rj = 
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Figure 16. Isosurfaces of dilatation, vorticity, temperature and concentration magnitudes in a 
128® subdomain covering (I 4277 )®. Yellow surfaces: 6 = —36'. (a) Green surfaces: uj — l.Soj'; (b) 
Green surfaces: Te = l.lTe'; (c) Blue surfaces: Co = —O.ACo', red surfaces: Co = l.QCo'. 




Figure 17. The p.d.f.s and conditional p.d.f.s of cosines of angles between scalar gradients and 
vorticity. —1 < 0 < 0: right triangles; —2 < 6 < —1: deltas; —00 < 6 < —2: squares; 0 < 0 < 1: 
left triangles; 1 < 9 <2-. gradients; 2 < 9 < + 00 : diamonds, (a) Temperature; (b) Goncentration. 
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Figure 18. Normalized dissipation spectra of kinetic energy and scalars. Inset: log-log plot of 

dissipation spectra. 


0.16, while that of concentration is about 2.35 and occurs at r/rj = 0.18. It indicates 
that the effectively smallest scale for concentration dissipation is around O.Ory, a bit 
smaller than the Kolmogorov scale. The position for the maximum dissipation spectrum 
of temperature is r/rj — 0.16, though the related OC scale is rjoc ~ lArj. So far the 
behind reason is unclear. In the inset we present the log-log plot. It shows that in the 
inertial range of O.I ^ r/rj < 0.18, the slope values of the dissipation spectra are about 
0.35, close to 1/3. In inertial range the spectra of kinetic energy and scalars follow the 
^-5/3 pQ.^^gj. therefore, the power law of the corresponding dissipation spectra should 
be 

In a stationary state, the kinetic energy and scalar dissipation rates oscillate in space. 
These oscillations are considered to cause held intermittency. It is therefore valuable to 
study the spatial structures of dissipation rates. Figure [T^ depicts the two-dimensional 
slices (x-y plane) through the dissipation helds simultaneously in the position of 2 ; = 7r/2. 
In order to observe the detailed structures at both small and large amplitudes, we plot 
the logarithms of dissipation rates, where the color scale is determined by the formulas 
as follows 

V'-x = logio (e/(e)), (5-1) 


= logio (e0/(e0)). (5-2) 

Here the colored range is —0.6 ~ 0.6. As dissipation increases, the color changes from 
blue to red. The kinetic energy dissipation rate is shown in Figure [T^a). We Hnd that the 
small-scale, high-dissipation regions, separated by the large-scale, low-dissipation regions, 
distribute discretely and randomly. These regions exhibit as sharp ridges. By contrast, 
the low-dissipation regions are more or less like deep canyons. Figure [T^b) shows the 
dissipation rate of temperature, it is consisted of the small-scale, cliff-like high-dissipation 
regions and the large-scale, ramp-like low-dissipation regions. Here the widths of the 
small-scale structures are even Hner, meaning less dissipation. In Figure 19 c), though 
the large-scale structures are similar to those shown in the former two panels, the small- 
scale structures are broader in widths, displaying as ribbons and thus indicating more 
dissipation. 

A common used method for quantifying intermittency is to compute intermittency 
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Figure 19. Two-dimensional slices of logarithms of the dissipation rates of kinetic energy and 
scalars at « = tt/2. (a) Velocity; (b) Temperature; (c) Concentration. 


parameter through the auto correlation of dissipation rate, namely 


(e<p(x)e^(x-hr)) ~ r 


(5.3) 


where fXip is the so-called intermittency parameter, and ip is for u, Te oi Co. The readers 
can find detailed definition in Wang et al. { 1999). In Figurej^we present the log-log plot 
of the auto correlations of e, ete and Cco, as functions of the normalized separation distance 
r/rj. There appear three linear regions for the auto correlations, in which /i„ = 0.14, 
Pte = 0.85 and p-co = 0.46. The relation of /i^ < Pco < ^J-te is in agreement with the 
fact that in our simulation, the intermittency of velocity, concentration and temperature 
increases sequentially. It should be noted that various methods (i.e. correlation method 
& variance method) may yield quite different values of intermittency parameter. 

Having analyzed the auto correlation, we shall now report on the correlation between 
the dissipation rates of kinetic energy and scalars. In Figure we plot the correlation 
coefficients between e and e^, as functions of the normalized separation distance r/ 77 , 
which are defined by 


Cie,e^) = 


(e- (e))(e 0 - (e,^)) 


((e-(e))^)i/2((e^_(e^))^)i/2 


(5.4) 


It shows that throughout scale ranges, both (e, ete) and (e, Cco) have positive correlations. 
The general behaviors of these correlation coefficients are that they decrease with scale 
and appear plateaus at the scales larger than Lf, where the level values are 0.05 and 
0.15, respectively. Note that the variation of correlation coefficient against scale in this 
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Figure 20. Auto correlations of the dissipation rates of kinetic energy and scalars, as functions 
of r/g, where the subscript ip denotes as u, te or co. Velocity: circles; temperature: squares; 
concentration: diamonds. 



Figure 21. Correlation coefficients between the dissipation rates of kinetic energy and scalars. 
Temperature: squares, concentration: diamonds. 


study is opposite to that in incompressible turbulence (Wang et al. 1999), in which the 
correlation coefficient grows with scale. In the range of r/rj ^ 2, C{e,ete) is larger than 
C(e, Cco), namely, the temperature is more closely connected with the velocity. Neverthe¬ 
less, as scale increases, this connection is weakened and is surpassed by the concentration. 
Furthermore, at small scales, C(e, Cco) is saturated at the level value of 0.25, whereas there 
is no saturation observed for C(e, Cte)- 

Finally, we turn to discuss the statistical dependence of scalar increment on dissipation 
rate. In Figure ]^ we plot the p.d.f.s of the normalized scalar increments conditioned on 
dissipation rates at the normalized separation distances of r/ry = 1, 3. For comparison, we 
also plot the unconditional p.d.f.s at the same separation distances. In the left panel, at 
each scale, the p.d.f. of temperature increment conditioned on the kinetic energy dissipa¬ 
tion rate, P{5rTe\e)^ is wider than that conditioned on the temperature dissipation rate, 
P{5rTe\ete). Moreover, P{5rTe\e) and P{5rTe\ete) become narrower as r/ry decreases, 
implying that in the dissipative range, the effects of dissipations make the large fluctua- 
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Figure 22. The p.d.f.s of the normalized scalar increments conditioned on the dissipation rates 
of kinetic energy and scalars. The solid lines with solid symbols represent the condition of e. 
The dashed lines with open symbols in (a) and (b) represent the conditions of ete and Cco, 
respectively. The dot-dashed and dot-dot-dashed lines represent the p.d.f.s of the normalized 
scalar increments a.t r/g — 1 and 3. 


tions of temperature happen less frequently at smaller scales. For each scale, the tails of 
P{SrTe) are longest, indicating that the dissipations suppress the intermittency of tem¬ 
perature increment. In the right panel, at each scale, P{5rCo\e) and P{SrCo) basically 
overlap each other, and are wider than P{6rCo\eco)- Unlike the 6rTe case, P{6rCo\e) be¬ 
comes wider but P{6rCo\€co) becomes narrower as r/g falls. Based on the above results, a 
conclusion is given that in compressible turbulence, the active scalar increment depends 
both on the kinetic energy and scalar dissipations, while the passive scalar increment 
mainly relates to the scalar dissipation. 


6. Cascades of active and passive scalars 

In this section, we employ a ”coarse-graining” approach ( Aluie|20TT Aluie et aZ.|2012 


Aluie|[20T3 ) to study the cascades of active and passive scalars. We begin with the defi¬ 
nition of a classically filtered field a;(x) 


i^(x) = J (firGi{r)a{y: + r), 


( 6 . 1 ) 


where G^(r) = G(r/l)/f is the kernel, and G(r) is a window function. The density- 
weighted Favre filtered field is then defined by 


a;(x) 


Pl(^) 


( 6 . 2 ) 


By the large-scale continuity and temperature equations, it is straightforward to derive 
the temperature variance budget for large scales as follows 

I il-piTei) +V-3i = -Ili-^i-Ai-Di + (6.3) 


Here J;(x) is the spatial transport of large-scale temperature variance, n^(x) is the 
subgrid-scale (SGS) temperature flux to scales < 1, 4)^(x) is the large-scale pressure- 
dilatation, A^(x) and DjiiPj are the viscous and temperature dissipations acting on scales 
> I, respectively, and £|°°^(x) is the temperature ’’energy” depletion caused by cooling 
function. The details of these terms are written as 
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Figure 23. Various terms except forcing in the stationary filtered equations of scalar 
variances, as functions of l/g. (a) Temperature; (b) Concentration. 


1 -- - ^ - 7 K d - 

^pTe Uj+pTe{Teuj Teuj) Te , 

(6.4) 

n = -p{TeUj - TeUj)-^Te, 

(6.5) 

$ = (7 — CjpOTe, 

( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 


Henceforth, we shall take the liberty of dropping subscript I whenever there is no risk of 
ambiguity. Similarly, the concentration variance budget for large scales is read as 

+V-31 =-Hi-Di + ef (6.9) 

where is the concentration ’’energy” injection from external forcing. The expressions 
of the spatial transport of large-scale concentration variance, J ^(x); the SGS concentra¬ 
tion flux to scales < I, n^(x); and the molecular dissipation acting on scales > I, Di{k), 
are 

Jj = -pCo Uj +pCo(Couj — Couj) — — Co , (6.10) 

—~p{CoUj — CoUj)— —Co, (6-11) 

OXj 

D=fP-Cco)\ ( 6 . 12 ) 

In Figure!^ we plot all the terms in Equations (6.3) and (6.9) except the forcing ones, 
as functions of the normalized scale l/r}. In the statistically homogeneous turbulence, 
the spatial transport of large-scale temperature variance vanishes, and thus, the value 
of the corresponding term shown in the left panel is basically zero. Throughout scale 
ranges, the SGS temperature flux and temperature dissipation are positive, causing the 
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Figure 24. Pressure-dilatation and its positive and negative components in the stationary 
filtered equation of temperature variance, as functions of l/rj. Total: squares; positive: deltas; 
negative: gradients. 



Figure 25. Normalized SGS fluxes of scalars and their components as functions of l/rj. (a) 

Temperature; (b) Concentration. 


temperature variance to decay. Contrarily, as scale increases, the magnitude of viscous 
dissipation declines quickly and approaches zero at the scales of Z/ry ^ 128, while that of 
pressure-dilatation first increases and undergoes a flat region in the range of 5 ^ Z/ry ^ 20, 
then decreases and reaches zero at large scales. Note that both the viscous dissipation 
and pressure-dilatation lead to the amplification of temperature variance. The above 
observation shows that the action of viscous dissipation is restricted at small scales, 
while the pressure-dilatation mainly takes place at the moderately large scales satisfying 
r] < I Lf. In the inset we enlarge the temperature dissipation and SGS temperature 
flux. It is found that the former limits itself at small scales and falls quickly as scale 
increases. By contrast, for the latter, there appears a plateau spanning the range of 
5 ^ Z/?y ^ 20, which demonstrates the existence of an inertial range for the temperature 
cascade. Throughout scale ranges in the right panel, the spatial transport of large-scale 
concentration variance vanishes as well. The molecular dissipation depletes concentration 
variance and occurs mainly at small scales. As scale increases, the SGS concentration 
flux first increases and undergoes a flat region of 30 ^ Z/ry < 80, then decreases at large 
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Figure 26. Isosurfaces of the compressive component of SGS scalar fluxes in a 128® subdomain 
(the filter width I = 5g), where the yellow and green surfaces are for = 0.25ete and 
= —0.25eco, respectively. 


scales. It implies that there is only direct cascade in the concentration transport, which 
is different from that in ID compressible turbulence ( Ni fc Chen||2012 |. 

The pressure-dilatation = (7 — ^)pfilTei does not contain any modes at scales 
< I, and thus, vanishes in the absence of SGS fluctuations. It only contributes to the 
conversion between the large-scale kinetic and internal energy, namely, if < 0, the 
energy is transported from kinetic to internal energy; if 0^ > 0, the process reverses 
( Aluie||2011 ). The picture of the negligible at small scales does not contradict to the 
fact of the rarefaction and compression motions appearing at all scales, which are the 
main property of compressible turbulence. Our result reveals that the high pressure- 
dilatation generated in the vicinity of small-scale shocklets will vanish after taking global 


averages, because of the cancelations between rarefaction and compression regions (Aluie 


et al. 20121. In Figure 24 we plot the pressure-dilatation and its positive and negative 
components, as functions of the normalized scale l/p. Obviously, both (4>^) and {^J) 
are high at small scales, however, when adding together, they basically cancel each other 
and make the outcome be small. 

In Figure!^ we plot the normalized SGS scalar fluxes and their components, as func¬ 
tions of the normalized scale I/t], In the left panel, there is a flat region in the range 
oi 8 ^ l/rj ^ 25, indicating the conservation of temperature cascade. The compressive 
component is significantly larger than the solenoidal component. This reveals that the 
cascade of temperature is mainly determined by the motions of rarefaction and compres¬ 
sion. By contrast, in the right panel, the plateau of the SGS concentration flux appears 
in the range of 35 ^ l/rj ^ 85, and the level value is close to unity. The magnitude of 
the solenoidal component is much higher than that of the compressive component. This 
means that the cascade of concentration is governed by the stretching and shearing of 
vortices. Furthermore, it is found that the compressibility of turbulence leads the com¬ 
pressive SGS concentration flux to be negative, and thus, trasfer from small to large 
scales. In a word, our results indicate that the cascades of the active and passive scalars 
are separately dominated by the compressive and solenoidal components of velocity. 
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Figure 27. The rescaled p.d.f.s of SGS temperature flux and its two components, (a) 

(b) P[inr/{eteii)]{i/v)-^-, (c) P[mrvietevm/g)-^ 


In Figurej^ the 3D isosurfaces of the compressive component of the SGS scalar fluxes, 
= 0.25ete and 11“° = —0.25eco, are displayed in a same 128^ subdomain, where the 
filter width is I = 5ry. It shows that the sheetlike surfaces of 11^°° are converged together, 
similar to the surfaces of dilatation shown in Figure [T6| This confirms the dominance of 
the rarefaction and compression motions in the temperature cascade. By contrast, the 
surfaces of 11“° are broken down into a large number of sheetlike fragments. 

The rescaled p.d.f.s of the SGS scalar fluxes and their components at I = 247], 32?7, ASrj 
and 64?7 are plotted in Figure 27 (temperature) and Figure (concentration). It shows 
that all the rescaled p.d.f.s of lT^°, 11^°® and 11^°° exhibit skewness toward positive side, 
therefore, the transfer of temperature flux is completely from large to small scales. By 
contrast, the rescaled p.d.f.s of 11“ and 11“® are positively skewed, whereas that of 11“° 
is negatively skewed. It indicates that the compressive component of the concentration 
flux transfers upscale. Here we point out that the p.d.f. tails of the SGS scalar fluxes are 
mainly contributed from the small-scale shocklets. Furthermore, we observe that in both 
figures, the rescaled p.d.f.s collapse to the same distribution for all I in the inertial range 
of 24 ^ l/rj ^ 64. By the multifractal theory (Benzi et al. 2008), the scaling exponents 
of the statistical moments of fl4^® and ffl“ as well as their components should saturate 
at high order numbers, where the values of exponents are equal to unity. This prediction 
reveals the statistical scale-invariant property of the active and passive scalar fluxes in the 
inertial range. Finally, the geometrical similarities of the rescaled p.d.f.s between (H^®, 
n^®°) and (n“, n“®), once again confirm the fact that the cascades of the active and 
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Figure 28. The rescaled p.d.f.s of SGS concentration flux and its two components, (a) 
p[mr/{t.ori)](ih)-^-, (b) P[mr (c) PiinrKe^ogWM)-^ 


passive scalars are respectively dominated by the compressive and solenoidal components 
of velocity. 


7. Summary and conclusions 

In this paper, a systematic investigation on the statistics of active and passive scalars 
in an isotropic compressible turbulent velocity was performed. The simulation was solved 
numerically using a hybrid method of a seventh-order WENO scheme for shock region and 
an eighth-order CCFD scheme for smooth region outside shock. The large-scale velocity 
and passive scalar forcings were added to drive the simulated system, and maintain 
it in a stationary sate, where the turbulent Mach number was Mf = 1.02, and the 
Taylor microscale Reynolds number was Re\ = 178. To explore the statistical differences 
between active and passive scalars, a variety of statistics including spectrum, structure 
function, probability distribution function, flow structure and cascade were investigated. 
Moreover, the effects of solenoidal and compressive processes of velocity on the scalar 
transports were described, of which some affect both active and passive scalars, while 
the others are shown to associate with only one scalar. 

The spectra of kinetic energy and scalars defer to the power law. The Kol¬ 

mogorov constant was found to be Ck = 2.06, while the OC constant for passive scalar 
is Cqq = 0.87, which is consistent with the typical values observed in incompressible 
turbulence. Compared to passive scalar, the transfer flux spectrum of active scalar shifts 
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towards higher wavenumbers, since the OC scale of temperature is larger than that 
of concentration. The local scaling exponents computed from the second-order structure 
function display flat regions for velocity and active scalar. Nevertheless, for passive scalar, 
it takes first a minimum of 0.61 and then a maximum of 0.73. Our results also confirmed 
that the velocity and passive scalar satisfy the Kolmogorov’s 4/5 and Yaglom’s 4/3-laws, 
respectively. However, there is no similar law found for active scalar. 

We then discussed the probability distribution function. First, the p.d.f.s of velocity 
and active scalar fluctuations are very close to Gaussian at small amplitudes and become 
super-Gaussian at large amplitudes, while that of passive scalar fluctuations arises slight 
oscillations at small amplitudes and changes into sub-Gaussian at large amplitudes. Sec¬ 
ond, the p.d.f. of velocity increment is concave and strongly intermittent at small scales. 
As scale increases, it approaches Gaussian. Further, by applying the Helmholtz decom¬ 
position to velocity, we found that at small scales, the p.d.f. for the solenoidal component 
of velocity is similar to that for velocity increment in incompressible turbulence, while 
the p.d.f. for the compressive component of velocity can be analogous to that for velocity 
increment in Burgers turbulence. As for the scalar increments, the p.d.f.s are symmetric 
and close to Gaussian at large scales. The shape for the p.d.f. of active scalar increment 
is concave, and that for the p.d.f. of passive scalar increment is convex, the same to 
that in incompressible turbulence. Third, the major contributions to the p.d.f. tails of 
flow gradients are caused by small-scale shocklets rather than large-scale shock waves, 
which make the values of the power-law exponents be —3.8 for velocity, —3.3 for active 
scalar, and —3.5 for passive scalar. Furthermore, our results showed that the skewness 
for velocity and passive scalar increments are positive, and approach zero as scale in¬ 
creases. Throughout scale ranges, in magnitude the skewness of passive scalar increment 
is always smaller than that of velocity increment, indicating a weaker degree of Gaussian 
departure. 

In terms of high-order statistics, the scaling exponent computed from structure func¬ 
tion gives the following relations: Zu,p > Zco,p > Zte,p, Zu,p > Zte,p > Ztem,p, and 
Zu,p > z,com,p > Zco,p-, where the last one is similar to that observed in incompressible 
turbulence. In addition, the computation based on the Helmholtz decomposition revealed 
that Zmtec,p is nearly identical to the mixed temperature-velocity scaling exponent shown 
in previous ID compressible turbulence simulation. 

We further investigated the field structures of scalars. It seemed that the active scalar 
has the ”ramp-cliff’ structures, which is usually regarded as production by stretching 
and shearing of vortices. By contrast, the passive scalar is dominated by the large-scale 
rarefaction and compression caused by shock fronts. The visualization of the isosurfaces of 
dilatation and scalars provides conclusions in two respects: (1) the sheetlike active scalar 
surfaces are wrinkled, distribute around shock wave and intersect with the dilatation 
surfaces in certain positions; and (2) the passive scalar surfaces are also sheetlike, and 
the opposite-sign value basically appears on each side of the shock fronts. The statistics 
of angles between vorticity and scalar gradients show that the scalar gradients prefer to 
be perpendicular to the vorticity. Compared to active scalar, the p.d.f. of passive scalar is 
much steeper, showing that the passive scalar is more frequently tangent to the vorticity. 

The analysis of dissipation field structures started from the dissipation spectra of 
scalars. We found that the maximum dissipation spectrum of passive scalar locates at 
higher wavenumbers than that of active scalar, given that the effectively smallest scale 
for passive scalar is about O.Qry, smaller than the Kolmogorov scale. In the range of 
O.I ^ r /?7 ^ 0.18, the dissipation spectra approximately follow a scaling, which can 
be derived from the Komlogorov theory. In order to understand the detailed structures 
at both small and large amplitudes, we computed the logarithms of the scalar dissipation 
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rates: tpte and 'tpco- It showed that ipte is composed of the small-scale cliff-like structures of 
high dissipation and the large-scale ramp-like structures of low dissipation. By contrast, 
the small-scale structures in ipco are broader in widths, and display as ribbons. The expo¬ 
nent parameter computed from the auto correlation of dissipation rate is a common used 
method for quantifying intermittency. Our results gave that Ute = 0.85 and fico = 0.46, 
meaning that the active scalar is much stronger in intermittency. As for the correlation 
coefficients between the dissipation rates of kinetic energy and scalars, it showed that 
C(e, Cte) and C(e, Cco) decrease with scale and emerges plateaus at scales larger than Lf. 
At small scales, C(e,eco) saturates at the level value of 0.25. The analysis was completed 
by the statistical dependence of scalar increment on dissipation rate. We found that the 
active scalar increment suffers influence from both the kinetic energy and active scalar 
dissipations, while the passive scalar increment mainly connects with its own dissipation. 

By employing a ” coarse-graining” approach, the scalar variance budgets for large scales 
were obtained. First, we observed that throughout scale ranges, the active scalar variance 
is increased by the large-scale pressure-dilatation and small-scale viscous dissipation, 
but is decreased by the SGS flux and temperature dissipation. The decomposition on 
pressure-dilatation showed that it has high values in the vicinity of small-scale shocklets. 
However, because of the cancelations between rarefaction and compression regions, the 
contribution vanishes after space average. As for passive scalar variance, it is depleted 
by the molecular dissipation at small scales, and is transported by the SGS flux from 
large to small scales. Second, the SGS active scalar flux, dominated by its compressive 
component, appears a plateau in the range oi 8 ^ l/rj ^ 25, implying the conservation of 
active scalar cascade. In spite of the compressive component is negative, the much larger 
positive solenoidal component leads the SGS passive scalar flux to transfer downscale as 
well, and defines the conservative cascade of passive scalar in the range of 35 < Z/ry ^ 85. 
As a result, in our simulation, the transport of passive scalar is mainly determined by 
the solenoidal component of velocity, which is similar to that in incompressible turbulent 
flows. Furthermore, the visualization for the isosurfaces of compressive-component SGS 
scalar flux at the filter width l/rj = 5 confirms the dominance of the motions of rarefac¬ 
tion and compression in the active scalar cascade. Third, it was found that in the inertial 
range, the rescaled scalar flux p.d.f.s collapse to the same distribution. According to the 
mutifractal theory, we predicted that the scaling exponents for the moments of and 
should saturate at high order numbers. This indicates that there are scale-invariant 
features for the statistics of active and passive scalars. In addition, the geometrical sim¬ 
ilarities of the rescaled p.d.f.s once again prove the aforementioned conclusion that the 
cascades of active and passive scalars are separately dominated by the compressive and 
solenoidal components of velocity. 

The current investigation reveals a variety of unique small-scale features of active 
and passive scalars in compressible turbulence, relative to the statistical properties of 
passive scalar in incompressible turbulence. Here we limit our study to the large-scale 
solenoidal forcing, a complementary numerical simulation in which both the solenoidal 
and compressive components of velocity are stirred will be performed ( Ni fc Ghen]20I5 ). 
Our future work will also address the effects of basic parameters such as the Mach and 
Schmidt numbers on the scalar transport in compressible turbulence (Ni]|20I5a|b|. 
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